HIGHER-RANK NUMERICAL RANGES AND 
COMPRESSION PROBLEMS 



MAN-DUEN CH0I\ DAVID W. KRIBS^'^ AND KAROL ZYCZKOWSKI^'^.e 

Abstract. We consider higher-rank versions of the standard nu- 
merical range for matrices. A central motivation for this investiga- 
tion comes from quantum error correction. We develop the basic 
structure theory for the higher-rank numerical ranges, and give 
a complete description in the Hermitian case. We also consider 
associated projection compression problems. 



1. Introduction 

In this paper we initiate the study of higher-rank versions of the 
standard numerical range for matrices. A primary motivation for us 
arises through the basic problem of error correction in quantum com- 
puting. Specifically, the development of theoretical and ultimately ex- 
perimental techniques to overcome the errors associated with quantum 
operations is central to continued advances in quantum computing. As 
it turns out, the numerical ranges Afc(T), for A; > 1, defined below are 
intimately related to this problem of "quantum error correction". In 
the paper fT] we give applications of the results from the present paper 
to this problem. 

Let T be an X matrix with complex entries. For A; > 1, define 
the rank-k numerical range of T as the subset Afc(T) of the complex 
plane given by 

(1) Afe(T) = {A G C : FTP = \P for some rank— A; projection Pj, 

where we use the term "projection" to mean "orthogonal projection". 
Observe that the numerical range of T is obtained as 

(2) Ai(T) = W{T) = {{T^ : G C^, || |^) || = 1}. 

In our analysis, it is desirable to explicitly find the scalars A and the 
associated projections P in Eq. (P). Thus, this "compression problem" 
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will be the focus of this paper. A search of the substantive literature 
on numerical ranges reveals connections with two lines of investigation. 
The "fc-numerical range" introduced by Halmos in [2J is the set of all A 
that satisfy A = Tr(PTP) for some rank-fc projection P. Evidently this 
set includes the set kAk{T), but the reverse inclusion does not hold in 
general. The "fcth matrix numerical range" studied by several authors 
consists of the set W{k : T) of all matrices X*TX, where X is an x A; 
matrix such that X*X = I. The higher-rank numerical ranges Afc(T) 
can alternatively be formulated as Afc(T) = W{k : T) fl {Xlk '■ A G C}. 
See 1^ 4i in] as examples of other entrance points into the literature 
on generalized notions of the numerical range. 

The rest of the paper is organized as follows. In the next section 
(§ 121) we discuss the basic structure theory for the sets Afc(T). In 
particular, we derive an explicit characterization of these sets for all 
Hermitian matrices. We state a conjecture and an open problem in 
the case of normal matrices. We discuss some lower dimensional cases 
in § 01 and in the penultimate section (§0)), we present a method for 
constructing the associated compression projections that captures all 
possible projections in the Hermitian case. In the context of quantum 
error correction, projections that correspond to elements of Afc(T), for 
A; > 1, must be explicitly identified. For instance, in the rank-two 
case, such projections correspond to quantum bits of information, or 
"qubits", that can be corrected after particular quantum operations 
act (see § Ell- 



in this section we investigate the basic structure theory of the sets 
Afc(T). We shall refer to elements of Afc(T) as "compression- values" 
for T, since A G Ak{T) if and only if the k x k scalar matrix Xlk is the 
compression of T to a /c-dimensional subspace. This means that T is 
unitarily equivalent to a 2 x 2 block matrix of the form 



where A is a k x [N — k) matrix, B is an {N — k) x k matrix, and C 
is an {N — k) x [N — k) matrix. Equivalently, T is a "dilation" of the 
scalar matrix A/fc, or, T — XI maps a fc-dimensional subspace into its 
orthogonal complement. 

The following set inclusions may be readily verified: 



2. Compression- Values 



(3) 




(4) W{T) = Ai(T) D A2(T) D . . . D A^{T). 

The following properties are also easily checked: 
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(z) Afc(aT + pi) = aAkjT) +P Va, /3 G C. 
(n) Afc(T*)=Afc(T). 
(m) Afc(T) C Afc(ReT) +iAfc(Imr). 
(zt;) Afc(T©5) D Afc(T)UA,(5). 

(v) Ak,+k,{T(BS)^Ak,{T)nAk,{S). 

The numerical range W{T) = Ai(T) is a non-empty, compact and 
convex subset of the plane that includes the spectrum of T. If T is 
normal, then W{T) is the convex hull of the eigenvalues for T. In 
particular, if T is Hermitian, then W{T) is the closed interval of the 
real line determined by the minimal and maximal eigenvalues of T. 
The higher-rank numerical ranges can, of course, be empty. But com- 
pactness still holds in general. The proof of the following result is 
elementary, hence we leave it to the interested reader. 

Proposition 2.1. Let T be an N x N matrix and letk>l. Then the 
rank-k numerical range Ak{T) forms a compact set. 

Now we give a description of the higher-rank numerical range for 
large values of k relative to N. 

Proposition 2.2. Let T be an NxN matrix and suppose that 2k > N. 
Then the rank-k nvmerical range Afe(T) is an empty set or a singleton 
set. If Ak{T) = {Ao} is a singleton set with 2k > N, then Aq is an 
eigenvalue of geometric multiplicity at least 2k — N. In particular, 
An{T) is non-empty if and only if T is a scalar matrix. 

Proof. Given 2k > N, assume that Afe(T) is non-empty, and contains 
Ao 7^ Ai. Let Pi be the corresponding rank k projections. Then 
the projection P = PqA Pi onto the intersection of the ranges of these 
two projections is non-zero and satisfies AoP = PTP = XiP. This 
contradiction shows that Afe(T) is a singleton set when it is non-empty. 
For the second claim, the equality P{T — XoI)P = imphes 

(5) T - Ao/ = (/ - P)(T - Aq/) + P(T - Ao/)(/ - P). 
Hence, rank(T - Aq/) < 2 rank(/ - P) = 2N - 2k, and so, 

(6) ker(T - Aq/) > N - {2N - 2k) ^2k-N. ■ 

In the normal case the previous result yields more detailed informa- 
tion for large values of k. 

Corollary 2.3. Let T be an N x N normal matrix and suppose that 
2k > N. Then the rank-k numerical range Ak{T) is an empty set or 
a singleton set. In fact, the case Afe(T) = {Aq} occurs if and only if 
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there is a {2N — 2k) x (2A^ — 2k) matrix Tq such that T is unitarily 
equivalent to Xohk-N © Tq, and Aq belongs to Ajv-fe(To). 

We now derive a general description of the rank-A; numerical range 
in the Hermitian case for arbitrary k. 

Theorem 2.4. Let A be an N x N Hermitian matrix with eigenvalues 
(counting multiplicities) given by ai < 02 . . . < oat and let k > 1 be a 
fixed integer with 1 < k < N . Then the rank-k numerical range Afe(>l) 

coincides with [ak,aN^k+i\ which is: 

[i) a nan- degenerate closed interval if < a^^k+i ? 
(u) a singleton set if ak = aN-k+i, 
(Hi) an empty set if ak > ojv-fc+i . 
Moreover, Ak{A) coincides with the intersection of the numerical ranges 
W{V*AV), where V runs through all isometrics V : C-^"*^"*"^ . 

Proof. Let A G Ak{A) and let be a rank-A; projection with PkAP^ = 
XPk- If V : is an isometry, then the subspace PkC^ 

and the range space VV*{C^) have non-zero intersection. Thus, there 
exists a unit vector e such that = Pk\i)) = VV*\ip). Let 
IV'') be the unit vector in C^~'^+^ given by \'4>') —V*\'^). Then we have 

(7) {V*AVil)'\'il)') = WIV') 

(8) = (PfeAPfeV'lV') = A(PfeVlV') = A. 

Hence we have shown that A belongs to W{V*AV). As V : C^-'=+^ 

was an arbitrary isometry, it follows that Ak{A) is contained in the 
intersection of all such numerical ranges 

Next, let : 1 < i < iV — A; + 1} be a fixed orthonormal basis for 
(^N-k+i Yet be an orthonormal basis for of eigenvectors 

for A corresponding to the eigenvalues Oi, . . . , otv- Consider two linear 
isometrics 14,^2 : C^-^+^ ^ defined by Vi{\i)) = ,V2{\i)) = 

\iJN-i+i)- 

Then V\*AVi and V2AV2 are operators on that are diagonal 

with respect to the basis and we have l^(y]*74yi) = [oi, aN-k+i] 

and W{V2AV2) — [ak, ojv]. It follows that 

(9) Ak(A) cf^W(V*AV) C W(V*AVi)f]W{V2AV2) 

V 

(10) = [ak,aM-k+i\- 

We complete the proof by showing Ak{A) contains the set [ak, aN^k+i] 
when Gk < aisf-k+i- Suppose first that aN+i-k > dk (and so 2k < N). 
Fix A in the interval [ak, aN+i-k]- We shall directly construct a rank-A; 
projection Pk such that PkAPk — XPk- Consider the set of k pairs 
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{ak-\-i-j,aN-k+j}, 1 < J < k. As a notational convenience we shall 
write {bj,b'j} for the ordered pair {ak+i-j,a^_k+j}, and so bj > b'-. 
(The following construction may be easily modified for any joint parti- 
tion of the sets {oat, . . . , aN-k+i} and {ofc, . . . , Oi} into ordered pairs.) 
We may write A, up to unitary equivalence, as a direct sum 

(11) A={®3 Aj) © 5, 

where each Aj is a diagonal 2x2 matrix with spectrum {6^,6^}, and 
B is either vacuous, or is the diagonal matrix with diagonal entries 
{afc+i, . . . , aN~k}- As A satisfies, 

(12) A G K, ajv-fe+i] C [6^., 6;.] = W{A,) V 1 < j < fc, 
we may find angles 9j such that 

(13) A = bj cos^ + b'j sin^ V 1 < j < k. 
Now define an orthonormal set of k vectors by 

(14) \(j)j) = cos 9 j\i)N-k+j) + sm.9j\i)k-j+i) ^1 <j <k, 

and the rank-/c projection onto the subspace spanned by these vec- 
tors; 

l0lX0l| + 102X021 + ... + |0fc)(0A:|. 

It follows that PkAPk = XPk- Indeed, observe that for 1 < j < k we 
have 

(15) = cos9i{AijN\(pj) +sm9i{A^i\4>j) 

(16) = aj^ cos 9i{ijj]\f\(f)j) + ai sin 9i{-ipi\(f)j) 

(17) = bi cos 9 1 cos 9j6N,N-i+j + b'l sin 9i sin 9j6j^i 

(18) = A5,,i. 

Similarly, {A(f)i\(f)j) = X6ij for I < i, j < k. 

The remaining case is characterized by the constraint X := ak = 
a^-k+i- If, in addition, a^_k+2 > ^fc-i, then we may split the sets 
{ttN, • • • , aAr_fc+2} and {afc_2, . . . , cii} into pairs as above, and similarly 
define k — 1 vectors \4>i), ■ ■ ■ , |0fc-i)- As the final vector we can take 
|0fc) := \^k), and define Pk = Y.''j=i l0iX0il- If «iv-fc+2 = Ofc-i, but 
aN-k+3 > CLk-2, then we will use {ipk) and l-^fc-i) as two of the vectors. 
This process may be continued, if required, to account for degeneracies 
in the spectrum of A around the eigenvalue a^, and construct a rank-fc 
projection which yields A G Ak{A). The result now follows. ■ 

For each real number r G M, we write [r] for the smallest integer 
n satisfying n > r. From Theorem 12. 4[ we see that if k < [A^/2] 
(equivalently 2k — 1 < N), then Ak{A) is non-empty for each x 
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a) ^1 ^2 H ^4 Ij-v a, 

-I ! \ ! ► ' -L- 

N=4 : :^A2^ : n=6 : 

A, ■ ■< 

Figure 1. The numerical range Ak{A) for a non- 
degenerate Hermitian operator A of size a) iV = 4, and 
b) = 6, with spectrum {oj}. 



A, 



_A2 
A, 



Hermitian matrix A. The following is an analogous result for a 
general non-Hermitian matrix. (This result can also be derived from 
Theorem 3.3 of [5 J.) 

Corollary 2.5. Let T be an N x N complex matrix. Let k be a positive 
integer satisfying k < [A^/4] (equivalently Ak — 3 < N). Then Afc(T) 
is non-empty. 

Proof. Write T = A + iB with A = A* andB = B*. Let b = 62^-1 be 
the (2/c — l)th smallest eigenvalue of B. By Theorem 12. 41 b G A2A;-i(-B); 
and so there is a projection P of rank 2k — 1 such that P{B — 6/)P = 0. 
Consider the {2k — 1) x [2k — 1) Hermitian matrix Aq given by the 
restriction of the compression PAP to the range of P. It follows from 
another application of Theorem 12.41 that Ak{Ao) is a singleton set {a}, 
where a is the fcth smallest eigenvalue of Aq. Hence there exists a 
projection Q < P such that rankQ = k and QAQ = aQ. Thus, 
QTQ = QAQ + iQBQ = (a + ib)Q and Afc(T) is non-empty. ■ 

The construction of projections that is described in the proof of the 
previous theorem will be further fleshed out in subsequent sections. It 
is perhaps appropriate to emphasize the most important non-trivial 
case of this result. Specifically, when 2k < N and the spectrum of 
A is non-degenerate. Theorem 12.41 shows that the rank-/c numerical 
range is the interval Ak{A) = [ak,aN-k+i] — see Fig. [U which shows 
generalized numerical ranges for A^ = 4 and A^ = 6. Also note that 
as an immediate consequence of Theorem 12.41 it follows that the sets 
Ak{A) are convex for all A; > 1 and Hermitian A. 

We finish this section by discussing the case of normal matrices. First 
note that property {Hi) above and Theorem 12.41 give a crude contain- 
ment result for Afc(T) for arbitrary T. Indeed, Afc(T) is a subset of the 
rectangular region in the complex plane {a + ip : a & Ajt(Re(T)), P G 
Afc(Im(T))}. In general we can obtain a more refined containment in 
the normal case. 
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Theorem 2.6. Let T be an N x N normal matrix and let k > 1 be a 

fixed positive integer. Then 

(19) Afe(T) C nr(cor), 

where T runs through all {N + 1 — k) -point subsets (counting multiplic- 
ities) of the spectrum ofT. 

Proof. The relevant parts of the proof of Theorem 12.41 can be easily 
extended to the normal case to verify the inclusion of Eq. ()19|) . ■ 

Remark 2.7. Observe that Theorem 12.41 shows the converse inclusion 
of Eq. ()19|1 holds in the Hermitian case. We believe this inclusion holds 
more generally, at least in the normal case, and we plan to undertake 
this investigation elsewhere. 

Conjecture 2.8. If T is an A^x normal matrix, then Afc(T) coincides 
with the intersection of the convex hulls co F, where F is an (A^-|- 1 — k)- 
point subset (counting multiplicities) of the spectrum of T. 

Verification of this conjecture would, of course, automatically imply 
that Afc(T) is convex, whenever this set is non-empty and T is normal. 
We state the general case as an open problem. 

Problem 2.9. Is Afc(T) a convex set whenever it is non-empty? 

As a consequence of Theorem 12.61 and the proof of Theorem 12.41 the 
conjecture can be seen to hold for A^ = 1, 2, 3, 4 and all values of k in 
each of these cases. Indeed, Theorem 12 . 6 1 shows that Afc(T) is contained 
in the desired set, and the construction of projections in the proof of 
Theorem 12.41 mav be adapted for A^ < 4. In each of these cases, Ak{T) 
is either the empty set, a singleton set, or an interval, and hence can 
never have interior. See Fig. |21for an illustration of some of the non- 
interval cases for A^ = 4. (We note that the A^ = 4 unitary case is 
explicitly worked out in [Ij.) 

The first open case is that of A^ = 5 and k = 2. The cyclic 5-shift 
is a good test example. This is the unitary U : defined 
on an orthonormal basis . . . , 1^5)} by U\Q = |^j+i(mod5))- The 

spectrum of U is given by Zn = exp{i^}, for n = 0, 1,2,3,4. Thus, 
A2(f/) is a subset of the pentagon shaped region depicted in Fig. 3. 
The arguments of (6) =^ (a) in Theorem 12.41 mav be used to show that 
A2(f/) contains the border points of this region, and also contains the 
centre A = 0. The problem is to determine if the rest of the interior 
points are included. 
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Figure 2. Examples of the numerical range A^iT) for 
a normal operator T of size a,b) = 4, and c,d) = 6, 
with non-degenerate complex eigenvalues {zi}. The set 
A2 is contained in the subset depicted (and equal at least 
in the case = 4). 




Figure 3. Numerical ranges Ak{U) for the cyclic 5-shift, 
with spectrum consisting of the fifth roots of unity, 
&) k — 1 and b) A; = 2, A2(C/) is contained in this set. 
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3. Eigenvalue- Pairing Construction 

The method presented in the next section shows how all the higher- 
rank projections may be obtained through a generalization of the "eigen- 
value pairing" approach used in the proof of Theorem 12 .41 For illustra- 
tion purposes, in this section we further discuss the pairing approach 
in some lower dimensional cases. 

First let us recall the k = 1 motivation for what follows 

below. If a AT > ... > ai are the eigenvalues of A = A* as above, 
then the numerical range of A is given by Ai{A) = [ai,aN]- (Assume 
A is non-scalar, so this is truly an interval.) Let be a choice of 
eigenvector for each aj. We may write a given A G Ai(A) as a linear 
combination A = aicl + a2cl + . . .+aNC%, where c = (q) are real scalars 
belonging to the (A^ — l)-dimensional simplex (YliLi ^1 ~ l)- Then the 
typical unit vector \(j)) that satisfies A = {A(f)\(j)), and the corresponding 
rank one projection, are given by 

N 

(20) |0) = ^e^^^c,|^,-) and P = 

i=i 

Observe that there are A^ variables Cj, with two constraints, and this 
gives N — 2 free parameters. There is also an additional A^ free phases 
from the choices 6j, for 1 < j < N. (Although it is A^ — 1 free phases 
up to a global phase allowed in the definition of P.) 

For instance, if {N,k) = (2,1), then Ai{A) = [01,02] where aj is 
an eigenvalue for the eigenstate In this case, A G [01,02] may be 
written as A = Oi cos^ l3i+o,2 sin^ Pi, and the angle Pi may be computed 
via the equation 

(21) cos2/3i = A^. 

Oi — O2 

The corresponding projection Pi is obtained as a "coherent mixture" 
of both eigenstates; 

(22) |0i) = e''^ cosAlV'i) + e''^ sin/5i|z^2), Pi = 

Next consider the case {N,k) = (3,1). Let A belong to Ai{A) = 
[oi, 03]. In this case, A may be obtained via the equation 

(23) A = oi cos^ Pi + 02 sin^ Pi cos^ P2 + 03 sin^ Pi sin^ P2. 

In this case, P2 = P2{^, Pi) depends on both A and Pi, and hence there 
is a one parameter family of solutions determined by Pi. The projection 
is given by Pi = |0i)(0i| where 

10,) = e''' cos/3i|V^i) + e''' sin/?i cos/32|V'2) + e''' sin/?i sin/?2|^3), 
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and we have three free phases {^^i, 6^2, 6^3} (two phases up to a global 
phase). In the case that A = 02, we may also use the solution Eq. (j22j) 
to find a vector |0i3) as a mixture of \ipi) and and then mix it 

with \il)2) to obtain 

|0l) = COS/32|V^2) + Sin/32|013)- 

Let us turn now to higher-rank projections obtained from the eigen- 
value pairing approach in the case iV = 4. The case of interest when 
iV = 4 is {N,k) = (4,2). The challenge occurs when A2{A) = [02,03] 
is a true interval. If we are given A G A2{A), we can consider all pairs 
{ai,ai'} that contain A. Here there are two possibilities: 

(i) {04,02}, {03,01}, 

(a) {04,01}, {03,02}. 

Of course, in the case of arbitrary N, there will be many more possible 
pairings. Now we solve the (2, 1) problem for each of the pairs sepa- 
rately. For instance, in the case of (i), we solve for f3i and (32 in the 
equations, 

f A = oi cos^ Pi + 03 sin^ Pi 
y \ = 02 cos^ /52 + 04 sin^ /32, 

and so, 

(25) cos^ Pi = — and cos^ Pi 



(26) 



Oi — O3 O2 — O4 

We then define coherent combinations of eigenstates grouped in pairs, 

\(f)^) = e''^' cos pi\tpi) + e''^^ sin Pil^s) 
102) = e'^^cosP2\iJ2) + e'^^smp2\iJ4). 

Then write P2 = |0iX0i| + |02)(02|, and it follows that PAP = XP. 

As noted above, this problem is equivalent to finding a unitary matrix 
U such that the matrix A' = UAU* includes a 2 x 2 block given by 
the scalar matrix Al2- In the case A = diag (oi, 02, 03, 04), observe 
that one choice for such a unitary is given hj U = OD, where D = 
diag (e*^^, e*^^, e*^^, e*^") and O is the orthogonal matrix given by 



(27) O 




sm/3i 
cos 132 sin (32 

cos /3i 
— sin 132 cos (32 



4. Higher-Rank Projections 

In this section we consider the problem of finding all possible rank-A; 
projections P associated with a compression- value A G Ak(T); i.e., to 
solve for the rank-/c projections P such that PTP = XP. We shall focus 
on the Hermitian case A = A*. Recall that in the proof of Theorem 
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and the discussion of the previous section we exphcitly constructed cer- 
tain famihes of projections to show that particular values of A belonged 
to A};(A). However, what we would like is a method for constructing 
such projections that captures all possibilities. Unlike the standard 
eigenvalue and eigenspace problem, in the generic case of this compres- 
sion problem there will be infinitely many projections. Indeed, even 
in the typical case for the numerical range = A.i[A) this is the 

case. But it is possible, and in fact easy, to write down such a method 
for the k = 1 case. There are of course more complications for k >2. 

First let us note that, while the eigenvalue-pairing approach con- 
structs a diverse set of projections, it is not sufficient to capture all 
projections associated with values of Ak{A). Indeed, even consider the 
k = 1 case of a Hcrmitian matrix A with spectrum {0,1,2}. Here, 
Ai(y4) = [0,2]. Let i = 0,1,2, be unit eigenvectors for the cor- 
responding eigenvalues. The eigenvalue-pairing approach for A = 1 in 
this case yields the family of projections P — IV'XV'I) where 

(28) |^) = _L(e^^i|^^)+e^^^|^2)). 

But the set of all projections P = \^){'4>\ such that {At/jltp) = 1 is the 
larger set given by unit vectors of the form 

(29) IV) = e^'^colVo) + cilV'i) + e''^co|V'2). 

For an arbitrary Hermitian matrix A, the rank one projections P = 
associated with values A G Ai(A) may be computed in the fol- 
lowing manner. Let Oi < ... < ajv be the eigenvalues for A, and let 
iVi), 1 < i < N, he a. choice of corresponding eigenvectors. Suppose 
we have a unit vector l'^) = X^^li Cj|'0i)- Then a simple computation 
shows that 

N 

(30) A = (A^IV) if and only if A = ^ai\ci\^. 

i=l 

This constructive condition characterizes the rank one projections 
associated with elements of the numerical range. Notice that there are 
infinitely many possibilities for such projections whenever A is not an 
eigenvalue for A. There is a corresponding characterization for arbi- 
trary k, though it is not constructive for k >2. Instead, in what follows 
we present a constructive, algorithmic approach to find all higher-rank 
projections associated with compression-values of Ak{A) for k > 2. 

Let A G Ak{A). By using a translation, we may assume that A = 
0. Let P+, Pq, and P_ be the projections onto the eigenspaces of 
A for respectively, the positive eigenvalues, the eigenvalue zero, and 
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the negative eigenvalues. First, we consider the case when there is no 
degeneracy in the spectrum of A; that is, A = is not an eigenvalue 
of A. Next choose a /c-dimensional subspace V+ of P+C^. Note that 
this is possible by Theorem 12 .41 and our assumptions A = G Ak{A) = 
[ofc, OTv-fc+i], and the non-degeneracy of the spectrum. By the same 
reasoning, P_C^ is at least fc- dimensional, and hence we may choose 
an isometry U : V+ ^ P_C^. Now we define a fc-dimensional subspace 
of C^; 

(31) V = {v+ + Uv+ -.v^ eV+}. 

Next define a /c-dimensional subspace W = f{A)V where 



(32) m-, , 

Observe that f{A)Af{A) = — P^, and hence Vfi, f 2 G V we have 

(33) {Af{A)v^\f{A)v2) = {f{A)Af{A)v^\v2) 

(34) = {(P^ - P_)v,\v2) 

(35) = {P+Vi\P+V2) - {P-Vi\P-V2) = 0. 

It follows that Pyv is a rank-A; projection such that Pyi;APy[) = 0. 

Now we show that every rank-/c projection P such that PAP = 0, can 
be written in the form P = Py^ as above. Let P be such a projection, 
and let V be the fc-dimensional subspace V = f{A)~^yV. Then for all 
f G V we have 

(36) = {Af{A)v\f{A)v) = {f{A)Af{A)v\v) 

(37) = {{P+- P-)v\v) = \\P+v\\^ -WP-vW^. 

In particular, this implies that the map U{P+v) = P-V determines a 
well defined isometry U : V+ — >^ V_, where V+ = P+V and V_ = P_V. 
Thus, V+ and V_ are both /c-dimensional and V is of the form given in 
Eq. (jSU, and hence P = Pw as claimed. 

We have presented a constructive method to obtain projections asso- 
ciated with the compression- values of Ak{A), in the case that there are 
no degeneracies in the spectrum of A. We have also shown that every 
such projection arises in this manner. Let us summarize the method. 

(i) Choose a fc-dimensional subspace V+ of P+C^. 
(ii) Choose a linear isometry f/ : V+ — P_C^. 
(Hi) Define the fc-dimensional subspace V = {v+ + Uv+ : E V+}. 
[iv) Let W = f{A)V. Then dim W = A; and PwAPw = 0. 

If there are degeneracies in the spectrum of A, the above method may 
be adjusted by including part of the subspace PqC^ in the subspaces 
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V+ and V_ as follows: As above, we want to construct all fc-dimensional 
subspaces V of such that 

(38) ((P+ - P^)v\v) =0 Vf G V. 

This can be accomplished since G Ak{A) = [ak,aN-k+i], and so 
k < dimPoC^ + dimP-tC^. Consider all possible pairs of non-zero 
integers {ki, ^2) with ^1 + ^2 = k, ki < dimPoC^, and 

(39) k2 < min{dimP_C^, dimP+C^}. 

Choose Vo as any fci-dimensional subspace of PqC^ and choose V+ as 
any fc2-dimensional subspace of P+C^. Let X : V+ ^ Cl^o)"*" H PqC^ 
be any operator and let ?7 : V+ — P_(C^) be any isometry and define 

(40) V = Vo + {v + Xv + Uv:v eV+}. 

Then W = f{A)V is a A;-dimensional subspace of C'^ with the desired 
properties. 

5. Concluding Remark 

We conclude by briefly discussing the mathematical context of the 
work [Ij, which includes applications of the present work to quantum 
computing. Every quantum operation £^ on a given quantum system 
is determined operationally by a set of operators {Ai} that act on the 
Hilbert space for the system via the so-called operator-sum representa- 
tion S{-) = J2i Ai{-)A*. (See |7| for a brief introduction to some of the 
mathematical aspects of quantum computing.) In the context of quan- 
tum error correction, the Ai are often called "error operators" . It is the 
effects of such operators that must be mitigated for whenever there is 
a transfer of quantum information determined by S. There are numer- 
ous strategies that have been, and are being, developed for this type of 
error correction. We go into detail on this subject in jT], but here we 
indicate how the mathematical conditions that characterize correction 
in the fundamental protocol for quantum error correction connects with 
the higher-rank numerical ranges. In the "standard model" for quan- 
tum error correction jSl IH], codes are identified with subspaces of the 
system Hilbert space, and "correctability" of a given code subspace C 
in terms of an error model S is shown to be equivalent to the existence 
of scalars A = [Xij) such that 

(41) PcA*A,Pc = K,Pc Vz,j. 

Here Pc denotes the projection of the system space onto C. Thus, the 
problem of finding correctable codes for a given error model S = {Ai} 
is equivalent to finding the compression-values inside the higher-rank 
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numerical ranges Ak{A*Aj), and V/c > 1, along with the corre- 
sponding projections. As indicated in P^, this problem may be reduced 
to a system of such problems for Hermitian or normal operators. 
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